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ABSTRACT 

In this contribution, we reexamine the radial oxygen abundance gradient in the strongly barred 
spiral galaxy NGC 3359, for which, using an imaging spectrophotometric technique, Martin & Roy 
detected a break near the effective radius of the galaxy. We have new emission line flux measurements 
of HII regions in NGC 3359 from spectra obtained with the Subaru telescope to further investigate this 
claim. We find that there are small systematic variations in the line ratios determined from narrow- 
band imaging as compared to our spectroscopic measurements. We derive and apply a correction to 
the line ratios found by Martin & Roy and statistically examine the validity of the gradient break 
proposed for NGC 3359 using recently developed metallicity diagnostics. We find that, with a high 
degree of confidence, a model with a break fits the data significantly better than one without it. This 
suggests that the presence of a strong bar in spiral galaxies can generate measurable changes in the 
radial distribution of metals. 

Subject headings: galaxies: individual (NGC 3359) — galaxies: abundances — galaxies: evolution — 
galaxies: ISM — galaxies: spiral 



1. INTRODUCTION 

Heavy metals are formed in stars and dispersed into 
the interstellar medium via mass loss processes leading 
to chemical enrichment. Spectral analysis of HII regions 
within gala xies can provide a measure of the chemical 
abundance (jPagel et al.lll979l ). The chemical abundance 
reflects the star formation history and the radial distri- 
bution of the chemical abundances in galaxies provides 
a key diagnostic for understanding galaxy formati on and 
evolu tion ( Vila-Costas & Edmunds 1992; Zaritsk y et al.l 
I1994D . A negative 0/H gradient has been established 
in a large number of normal and barred spiral galax- 
ies (lSearlelll971l : iGarnett fc Shieldslll987HBresolin et al.l 
I2OO9I . ainong others). This is thought to be a result of 
greater gas processing in the center of these galaxies. 

It has been found that the 0/H abundance in barred 
galaxies is considerably more shallow than the gradient 
of normal spiral galaxies an d this flattening is correlated 
to the strength of the bar (jMartin fc Rovlll994f ). This 
has been attributed to ho mogenization of th e abundances 
by large scale gas flows (Pagel et al." '1979t lAUoin et al.l 
[1981). Simulations by Friedli et al. (19941) have demon- 
strated that within the corotation radius gas accretion 
leads to increased star formation, resulting in an in- 
crease in the chemical abundance, while large-scale gas 
motions leads to a mixing outside the corotation radius. 
This mixing homogenizes the outer part of the galaxy 
resulting in a break at the corotation radius in the radial 
abundance gradient. This break is t herefore taken as an 
indicator of a recently formed bar (|Rov fc WalshI 119971 : 
IFriedlil[T999t IConsidere et al.l[20nnl) . 

NGC 3359 is a barred spiral ga l axy cl assifled as a 
SBc(rs)II by Ide Vaucouleurs et al.l ()1991[ ). The bar 
ratio has been determined to be b/a = 0.3 , mak- 
ing N G C 3359 a very strongly barred galaxy (|Martinl 
[T99I . iMartin fc Rol (fl995l herafter MR95) have de- 
termined the abundance gradient for NGC 3359 us- 
ing the [Oil I]/ HP and [Nil] /[Oil I] abundance indi- 
cators obtained from monochromatic images of nebu- 



lar emission lines of 77 HII regions. MR95 measure a 
steep slope, determined from a linear regression, for the 
0/H gradient inside the corotation radius of 4.7 kpc 
{Alog{0/H)/AR = -0.070 ± 0.010 dex kpc^^) and a 
flattening outside {Alog{0/H)/AR = -0.006 ± 0.018 
dex kpc^^) which they attribute to the homogenization 
of the abundance due to gas mixing. To date there exists 
only a few examples in the literature of such breaks and 
this break in the slope observed by MR95 in NGC 3359 is 
one the clearest examples. Though models suggest that 
the formation and presence of bars result in a break in 
the abundance gradient, the cause and validity of the 
breaks are still not well established obs ervationally. 

It has been argued bv lPilvuginl ()2003f ) that these bends 
may be attributable to systematic errors or improper 
use of the relationship between abundances and strong 
line intensities. He asserts that abundance indicators 
constructed from the oxygen strong-lines reproduces the 
metallicities for the high-abundance HII regions, but 
overestimates the metallicity in the low-abundance pe- 
ripheral regions of galaxies leading to a false bend in the 
gradient. Moreover, he finds that a false bend in the 
abundance gradient may result from a bend in the ion- 
ization parameter gradient. 

We reexamine the abundance gradient of NGC 3359 
and apply recently developed diagnostics appropriately 
calibrated and suited for this study. In §2 we present 
our data and methods, in §3 we compare our spectro- 
scopic line ratio measurements to the photometric mea- 
surements of MR95. In §4 we compare a model with two 
lines and a break with a single line model by performing 
a likelihood ratio test. In §5 we examine the ionization 
parameter and in §6 we present a brief discussion and 
our conclusions. In this work, we use a galactic distance 
of 14.4 Mpc for NGC 3359 determined from the radial 
velocity of 1008 km s"-'^, using Hq = 70 km Mpc~^. 
FolowingHal (|1986[ ). we adopt 51° and 172° for the in- 
clination and position angle of NGC 3359, respectively. 
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Fig. 1.— Ha image of NGC 3359 taken with FOCAS on-board 
the 8.2m Subaru telescope. A prominent bar, showing a chain of 
HII regions, dominates the Ha emission. Two distinct spiral arms 
can be seen with numerous HII regions highlighting the spiral struc- 
ture. The 19 regions that have been observed spectroscopically and 
that meet our S/N cut are identified in the image. The image is 
oriented with North at the top and East to the left. 

2. DATA 

NGC 3359 was observed using the Faint Object Cam- 
era and Spectrograph (FOCAS) on-board the 8.2m Sub - 
aru telescope on Mauna Kea (jKashikawa et al.1 12002D . 
Multi-object spectroscopy of 26 HII regions ranging from 
the center out to the periphery of the galaxy was taken. 
We obtained the blue portion of the spectrum from two 
30-minute exposures using the 300R grism in second or- 
der, and the red portion of the spectrum from two 20- 
minute exposures using the VPH600 grism. The obser- 
vations were made under photometric conditions with a 
seeing of 1.7". The nominal spectral coverage for the 
blue portion is typically 3700 — 5100 A, although in a 
few spectra the [O//] doublet centered at 3727 A is ei- 
ther not observed due to restricted wavelength coverage 
or lost in the noise. The nominal spectral coverage for 
the red portion of the spectra is 6300 — 8200 A. The 
spectral resolution of these spectra is ~ 5 A. 

Figure [T] shows an iJa image of NGC 3359 obtained 
as part of the preimaging for this project. The intense 
star formation in NGC 3359 is traced out by its ^ 100 
HII regions. There is a prominent bar with a misaligned 



TABLE 1 
HII Region Properties 



ID 




R/R25 


R (kpc) 




E(B-V) 


1 


1.17 


0.47 


7.18 


2.70 


0.52 


2 


1.19 


0.48 


7.31 


2.43 


0.28 


3 


1.35 


0.54 


8.29 


2.50 


0.35 


4 


1.36 


0.55 


8.38 


2.41 


0.27 


5 


0.86 


0.34 


5.26 


2.58 


0.42 


6 


0.47 


0.19 


2.91 


2.77 


0.58 


7 


0.82 


0.33 


5.05 


2.75 


0.56 


8 


1.11 


0.45 


6.85 


2.62 


0.45 


9 


1.80 


0.72 


11.0 


2.82 


0.62 


10 


0.16 


0.06 


0.95 


2.47 


0.33 


11 


0.69 


0.28 


4.25 


2.43 


0.29 


12 


0.57 


0.23 


3.49 


2.41 


0.27 


13 


1.42 


0.57 


8.73 


2.48 


0.33 


14 


1.45 


0.58 


8.92 


2.33 


0.19 


15 


1.49 


0.60 


9.15 


2.12 


0.00 


16 


1.57 


0.63 


9.63 


2.45 


0.30 


17 


1.64 


0.66 


10.0 


2.20 


0.07 


18 


1.70 


0.68 


10.4 




0.33 


19 


2.39 


0.96 


14.7 




0.33 



Note. — Properties of the 19 HII regions from spectroscopic 
line flux measurements for NGC 3359. The first column gives the 
ID corresponding to the HII regions labeled in Figure[T] Column 2 
and 3 give the deprojected radial distances relative to the effective 
{Reff = 1-4 arcmin) and 25th B magnitude isophotal (-R25 = 3.6 ar- 
cmin) radii, respectively, taken from [de Vaucouleurs et al.l 1)19911 ). 
Column 4 gives the deprojected galactocentric distances in kpc as- 
suming a distance of 14.4 Mpc. Column 5 and 6 give the Balmer 
decrement determined from the H BJH^ ratio and and extinction 
inferred from the reddening law of lCardelli et aO I I1989I) . 

chain of HII regions along with two spiral arms. 

The spectroscopic data were reduced and ID spectra 
were obtained using a suite of IDL routines developed 
specifically for these data. For each region, we perform 
an individual sky level subtraction by averaging the off- 
source pixels in each slit. The data were flux calibrated 
in IRAF^ using spectra of four observed standard stars. 
The fluxes for the emission lines were determined by sum- 
ming all the pixels in the line bandpass while subtracting 
out the continuum determined from a linear fit to a re- 
gion immediately outside the bandpass. 

The errors in the emission line flux, cr;, 
were determined using th e equation given in 
IGonzalez-Delgado et all (|1994D : 

ai = a,iVi/2[l + EW/{N^)Y/^. (1) 

The error in the continuum, a^, is measured in the re- 
gion near the emission line over which the linear fit to 
the continuum was performed. N is the width of the re- 
gion in pixels over which the flux in the emission line is 
measured. EW is the equivalent width of the emission 
line and A is the spectral resolution in A pixel" ^. The 
first term in the equation accounts for the uncertainty 
in the continuum near each line by measuring the rms of 
the continuum and multiplying by the width of the 'base' 
of each emission line. The second term accounts for the 
Poisson error in the line intensity, which is expressed here 
through the equivalent width and is added in quadrature. 
For this work, where the S/N in our lines of interest is 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 2 

HII Region Emission Line Strengths 



ID 


[Oil] 


[OIII] 


[OIII] 


[OIII] 


[Nil] 


[SII] 


[SII] 


12+log(0/H) 


12+log(0/H) 




3727 


4363 


4959 


5007 


6583 


6717 


6731 


(N2) 


(03N2) 


1 


082 ± 3d 


■ ■ ■ 


60 it 4.2 


184 ± 24 


31 it 2.7 


36 it 2.4 


25 it 2.3 


o on 1 n r\o 

8.30 ± 0.02 


o o ^ In nn 

8.34 ± 0.02 


2 


o /I 1 fin 

245 lb 2y 


1.7 it 0.3 


1 on 1 c o 

130 ± 6.8 


ion 1 or\ 

389 ± 20 


18 it 1.1 


20 it 1.1 


16 it 1.7 


o on 1 n m 

8.20 ± 0.01 


o -\ c In m 

8.16 ± 0.01 


3 


205 it 31 


2.8 it 0.6 


138 it 7.4 


406 ± 20 


10 it 0.7 


9.3 it 0.5 


6.8 it 0.5 


Oil In no 

8.11 ± 0.02 


o no 1 n ni 

8.08 ± 0.01 


4 






44 it 2.9 


135 it 7.4 


36 it 2.6 


38 ± 2.7 


27 it 2.1 


o o o 1 n no 

8.33 ± 0.02 


o /I n 1 n ni 

8.40 ± 0.01 


5 


231 ih 18 




58 it 3.1 


173 ± 9.2 


43 it 2.3 


32 it 1.7 


24 it 1.5 


8.38 ± 0.02 


o QO _l_ n ni 
8.39 ± 0.01 


6 


314 it 63 




13 it 1.7 


41 it 3.4 


53 it 3.0 


38 it 2.3 


27 it 1.6 


8.44 ± 0.02 


o £jo _l_ n ni 
8.62 ± 0.01 


7 


2o7 ± 23 




71 it 3.7 


Z\Ji ±10 


26 it 1.4 


22 it 1.2 


16 it 0.8 


Q O'/ _l_ n ni 
8.2 / ± 0.01 


o Qn _l_ n ni 
8.30 ± 0.01 


8 


490 it 64 




29 it 2.2 


74 it 4.4 


36 it 2.4 


32 it 2.3 


23 it 1.7 


8.33 it 0.02 


8.49 it 0.01 


9 


410 it 59 




78 it 6.8 


223 it 17 


22 it 4.1 


19 it 4.8 


18 it 4.6 


8.23 it 0.04 


8.26 it 0.03 


10 


116 it 11 




17 it 2.5 


60 it 4.0 


71 it 3.9 


38 it 2.0 


28 it 1.6 


8.54 it 0.03 


8.61 ±0.01 


11 






51 it 2.6 


152 it 7.7 


39 it 2.2 


35 it 1.8 


25 it 1.3 


8.35 it 0.02 


8.40 it 0.01 


12 


272 it 15 




45 it 2.4 


135 it 7.0 


40 it 2.2 


27 it 1.5 


19 it 1.0 


8.36 it 0.02 


8.42 it 0.01 


13 


270 it 17 


2.7 it 0.4 


91 it 4.8 


270 it 13 


20 it 1.1 


27 it 1.4 


19 it 1.4 


8.22it0.01 


8.23 it 0.01 


14 


382 it 31 




64 it 3.6 


191 it 10 


32 it 2.0 


40 it 2.3 


28 it 1.7 


8.31 it 0.02 


8.34 it 0.01 


15 






111 it 6.3 


331 it 17 


23 it 7.1 


38 it 6.6 


33 it 8.6 


8.24 it 0.24 


8.22 it 0.05 


16 






85 it 4.5 


250 it 12 


28 it 1.6 


41 it 2.2 


31 it 1.7 


8.28 it 0.02 


8.28 it 0.01 


17 


332 it 67 




48 it 4.5 


156 it 9.3 


29 it 3.9 


46 it 4.5 


40 it 4.4 


8.29 it 0.03 


8.35 it 0.02 


18 






62 it 11 


203 it 20 


38 it 12 


96 ± 30 


79 ± 36 


8.35 it 0.18 


8.35it0.11 


19 






102 it 12 


287 it 22 


17 it 5.9 


21 it 5.3 


13 it 2.6 


8.19it0.39 


8.19it0.11 



Note. — Emission line strengths and metallicities of the 19 HII regions from spectroscopic line flux measurements for 
NGC 3359. Column 1 gives the ID corresponding to the HII regions labeled in Figure [T] Columns 2-8 give the line flux relative 
to Hp of 100. Column 9 a nd 10 give the oxygen abundance, 12+log(0/H), derived from the A'^2 and 03N2 line ratios calibrated 
bv lPettini fc Page!! pOOl V The errors in the oxygen abundance are determined from propagating the uncertainties in the line 
ratios using a Monte Carlo method and do not reflect the systematic uncertainties of the strong line methods for abundance 
determinations. IKewlev & EUisor] II2008I ) provide a more detailed discussion of the systematic uncertainties. 



quite high, we consider this adequate for characterizing 
the Hne flux uncertainties. We account for additional 
uncertainties from the flat fielding, flux calibration and 
scaling due to the spectra being obtained through differ- 
ent grisms by adding in quadrature relative uncertainties 
of 1%, 4% and 3%, respectively. 

The data were der eddened using the reddening law of 
iCardelli et all ()1989D with the E(B-V) value determined 
from the Balmer decrement using the H/3/H"/ ratio. For 
Case B recombination (T p = 10*^ K, Up = 100 cm ~^), 



this ratio is equal to 2.14 (jHummer fc Store"vlll987[ ). For 
regions 18 and 19 we were not able determine the Balmer 
decrement due to non-detection of the H'y line. For these 
regions we have adopted the median value of the extinc- 
tion. After dereddening the line fluxes, the spectra ob- 
tained in the two grisms were scaled to have the Ha/H/3 
ratio consistent with Case B recombination. We select 
our final sample of 19 HII regions by imposing a min- 
imum signal-to-noise (S/N) ratio cut of 3 for the line 
r atios of interest. 

'Pettini fc Pagell ((2004) derive an empirical method for 
estimating oxygen abundances based on the [Nil]/ Ha 
and {[OIII]/ H {[N II]/ Ha) line ratios. They param- 
eterize the relationship between \N II]/ Ha and 12-1- 
log{0 / H) with a third degree polynomial given as 

12 + log{0/H) = 9.37 + 2.03■iV2 + 1.26•7V22 + 0.32■iV2^ 

(2) 

where N2 — log{[N II]/ Ha). Additionally, they derive 
a linear relationship between {[OIII]/ H P) / {[Nil]/ Ha) 
and 12 + log{0/H) given by 



12 + log{0/H) = 8.73 - 0.32 • 0?>N2, 



(3) 



where 03iV2 = log{{[OIII]/ H p) / {[Nil]/ Ha)}. These 
diagnostics are calibrated using HII regions where metal- 
licity has been determined from either direct measure- 
ments of the electron temperature or by detailed pho- 



toionization modeling. The diagnostics using the N2 and 
03N2 line ratios are monotonically increasing with O /H 
and are not double valued like many of the diagnostics 
using strong oxygen lines (i.e. R23). Additionally, these 
diagnostics are calibrated over the full metallicity range 
observed in the HII regions of NGC 3359. The la dis- 
persion of the N2 and 03N2 calibrations are 0.18 and 
0.14 dex, respectively. 

Table [1] lists some of the properties of the HII regions 
identified in Figure [1] and Table [2] gives the emission 
line strengths and metallicities. The error estimates for 
metallicities given in Table 2 only reflect the statistical 
uncertainties, the systematic uncertainties in the abso- 
lute metallicities are considerably larger. However, we 
note that for the metallicity analysis performed in this 
study, we only req uire relative accuracy wh ich the diag- 
nostics do provide ^Kewlev fc Ellisonll2008[ ). 

3. SPECTROSCOPIC VS. PHOTOMETRIC LINE RATIOS 

In this section we compare our spectroscopic line ra- 
tio measurements to those determined from narrow-band 
photo metry given in Table 4 of MR95. iDutil fc Rovl 
(|2001f ) point out some of the hazards of using imag- 
ing spectrophotometry for studying chemical abundances 
across extended objects such as galaxies. They demon- 
strate that the main limit to their method is the sub- 
traction of the underlying continuum. We compare our 
high S/N spectroscopic measurements with robust mea- 
surements of the underlying continuum with the spec- 
trophotometric measurements of 1VIR95 in order to asses 
any systematic bias that may result from photometric 
determinations . 

We match HII regions in common with ]V[R95 by eye, 
comparing our image shown in Figure[T]to a similar figure 
(Figure 3) in MR95. We note that for 2 of the 17 matched 
HII regions there is ambiguity owing to the high density 
of labeled points in 1VIR95. This ambiguity is inconse- 



4 




Fig. 2. — The photometric measurements of MR95 compared with our spectroscopic measurements. The figure shows the [OIII]/Hf3 
(left), [NII]/[OIII] (center) and [NII]/Ha (right) line ratios. The data points are the 17 regions for which there is both a photometric and 
spectroscopic line ratio measurement. The solid line in each of the plots is the one-to-one agreement and the dashed line is a least-square fit 
to the data. The error bars are the statistical uncertainties taken from MR95 (x-direction) and this study (y-direction) and do not reflect 
any of the systematic uncertainties. 



quential for one of the regions because the hne ratios for 
the two possible matches are consistent to within the er- 
rors. For the second case, we select the region that is 
most consistent with the observed relationship between 
the spectroscopic and spectrophotometric measurements 
inferred from the other 16 HII regions. The selection 
has very little effect on the analysis and no effect on the 
conclusions. 

Figure [2] shows the comparison for the [OIII]/ H (3 
(left), [Nil]/ [Oil I] (center) and [Nil]/ Ha (right) line 
ratios. The errors on the MR95 data are taken from 
their measurements. The 17 HII regions span the full 
metallicity range covered by the MR95 data. The solid 
line in each of the plots is the one-to-one agreement and 
the dashed line is a least-square fit to the relationship. 
The dispersion around the fitted relation is 0.12, 0.18 
and 0.10 for [OIII]/HP, [NII]/[OIII] and [NII]/Ha, 
respectively. The [NII]/[OIII] line ratio has the largest 
dispersion which we attribute to the additional scatter 
resulting from uncertainty in the reddening correction 
d ue to those lines bei ng ~ 1600 A apart. 

iDutil fc Rovl (|2001| ) show that the spectrophotometric 
method is limited by the continuum subtraction. The 
systematic deviation of the MR95 data from our data, 
highlighted by the linear fits in Figure [21 are consistent 
with a systematic error in the subtraction of the contin- 
uum. We expect a systematic error in the continuum 
subtraction to most strongly effect the weaker lines. For 
the ratios shown in Figure [21 the weaker line is -ff/3 for 
the [OIII]/Hl3 line ratio (left panel) and [Nil] for the 
[Nil] /[Oil I] and [Nil] /Ha line ratio (center and right 
panel, respectively) . The deviation observed is consistent 
with an overestimate of the continuum leading to an un- 
derestimate of the line strengths in the MR95 data. 

We derive a linear correction for each of the line ratios 
of the photometric data of MR95 from the 17 measure- 
ments displayed in Figure [2] by minimizing the square of 
the residuals, weighted by the statistical errors, between 
the two determinations and apply this correction to all 
77 line ratio measurements made by MR95. Hereafter, 
when referring to the data of MR95, we mean to refer to 
these corrected data. 

The linear fits shown in each of the three panels of 



Figure [2] were performed using the routine fitexy.pro in 
IDL. The routine does a linear least-square approxima- 
tion using the errors in both the x and y data. The fits 
are given by: 

SosHp = (-0.08 ± 0.02) + (1.32 ± 0.08)Po3H/j, (4) 
Sn203 = (-0.01 ± 0.02) + (1.24 ± 0M)Pn2O3: (5) 
and 

Sn2 = (0.12 ± 0.02) + (1.21 ± 0.07)Pjv2- (6) 

Here S and P denote the spectroscopic and photometric 
data, respectively. The subscripts in each of the equa- 
tions denote the line ratio in consideration. The errors 
in the linear fits in equations [31 - [H refiect measurement 
uncertainties but do not account for systematic uncer- 
tainties. 

We note that, to within the scatter, the line ra- 
tios appear to have only linear systematic differences. 
Furt hermore, the ab undances determined from the 
Pettini fc Pagell (|2004D diagnostics are nearly linear in 
metallicity with respect to the line ratio used. There- 
fore, it is extremely unlikely that an abundance gradient 
break could spuriously result from the use of the photo- 
metric as opposed spectroscopic line ratios. 

We conclude that given the greater systematic uncer- 
tainties in the spectrophotometric method of line ratio 
measurement, our spectroscopic data provide a more ro- 
bust measure of the line ratios used in ionized gas di- 
agnostics as compared with the monochromatic imaging 
data of MR95. However, the strong correlation and rel- 
atively small dispersion of MR95 with our data give us 
confidence in the corrected data of MR95 and its use in 
estimating abundances. 

4. ABUNDANCE GRADIENT BREAK 

The oxygen abundance gradient in NGC 3359, first in- 
vestigated by MR95, offers one of the clearest examples 
of a break in the radial gradient found for barred spi- 
ral galaxies. In this section, we reanalyze the data of 
MR95, but employ a more robust statistical approach 
for establishing the break. The data of MR95 are consis- 
tent with our observations and we have recalculated the 
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Fig. 3.— The abundance gradient for NGC 3359 derived from \NII]/Ha (left) and [0///]/[Af//] (right). The black data points are the 
metallicities determined from the N2 and 03N2 diagnostics of Pcttini & Pagel (2004) using the line ratios of MR95. The diamonds are 
the metallicities determined for 19 regions from our spectroscopic measurements of the line ratios using the same metallicity diagnostics. 
The solid line in each panel is a linear model fit to the MR95 data and the dashed lines in each panel is a model with a break fitted to the 
MR95 data. 



m etallicities using t he em pirically calibrated diagnostics 
of |Pettim & Page! (2004"). The N2 ratio does not suf- 
fer from the degeneracy found in the oxygen line ratios 
due to its monotonic behavior with respect to metallic- 
ity. Additionally, due to the very close spacing of the the 
two lines involved, the ratio does not require any correc- 
tions for reddening. We also calculate metallicities based 
on the 03N2 line ratio, though this diagnostic yields a 
considerably larger scatter in 0/H. 

Figure [3] shows the abundance gradient derived from 
the ^^2 (left) and 037V2 (right) diagnostics. The black 
data points are the metallicities determined from the 
data of MR95. The diamonds are the metallicities for 
19 HII regions derived from our spectroscopic observa- 
tions. The solid line in each plot is a least square fit to 
the data from MR95 while the dashed lines are a linear 
model with a break fit to the same data. The model with 
a break has five parameters, four of which define slopes 
and intercepts of two lines, while the fifth parameter de- 
fines the location of the break in units of the effective 
radius. MR95 take the break to occur at an effective 
radius of R^ff = 1. 

In order to establish the significance of a model with 
a break in the gradient we perform a likelihood ratio 
F-test. This test is appropriate when the distribution of 
the statistic is a F-distribution under the null hypothesis. 
This test is most commonly used when comparing nested 
statistical models such that the model with fewer param- 
eters, model 1, can be reproduced by fixing parameters 
in the model with greater number of parameters, model 
2. In our case, model 1 (linear model) can be reproduced 
by setting the break to be at Reff > Rmax, where Rmax 
is the maximum effective radius of the data, making the 
second line inconsequential when determining the sum of 
the residuals squared in model 2 (model with a break). 



Model 2 will always fit as well, if not better than model 1. 
The purpose of performing the F-test is to judge whether 
the decrease in the sum of the residuals squared is sig- 
nificant given the increase in the number of parameters. 

One can calculate the F-statistic for the null hypothesis 
that model 2 does not give a significantly better fit than 
model 1. The F-statistic is given by 



N-P2 

Here, 

^.^^0^^^ (8) 

where Xi — rrii are the residuals between the model and 
the observations and af are the errors in the observa- 
tions. P is the number of parameters and the subscripts 
denote the model. N is the number of data points. Our 
derived F values are 8.5 and 13.6 for the N2 and 03A^2 
diagnostics, respectively. Under the null hypothesis, F 
has a F-distribution with (P2 — Pi, TV — P2) degrees of 
freedom. We can derive the probability at which we can 
reject the null hypothesis by comparing our calculated F 
value with a critical value at some false-rejection proba- 
bility. For F distributed with (3, 72) degrees of freedom, 
the critical value is 8.1 for a false-rejection probability of 
0.0001. The calculated F values suggest that the prob- 
ability that we would falsely reject the null hypothesis 
is less than 0.01%. To put it in terms of a confidence, 
we can say with over 99.99% confidence that the null 
hypothesis is rejected, meaning that model 2 is a signifi- 
cantly better than model 1 at fitting the data. A F-test 
performed on our spectroscopic data alone gives a confi- 
dence of 75% and 65% for the A^2 and 037V2 diagnostics. 
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using the calibration of iKobulnickv fc KewlevI ()2004l) . 
The hne ratios are defined such that 



R23 = log{{[OII] + [OIII])/H/3} (9) 



and 



2.0 



Fig. 4. — The abundance (top panel) and ionization parameter, 
q, (bottom panel) derived fro m the R23 and 032 lin e ratios, re- 
spectively, using the scheme of lKobulnickv KewlevI (2004). The 
ionization parameter is in units of cm s~^. The upper/lower 
branch degenerac y is b roken using the [N I I]/[OI I] line ratio 
pewlev fc £11180312003 ). 

respectively. 

From visual inspection of the left panel of Figure |3] it 
can be seen that the gradient does not appear to be con- 
tinuous. The greater scatter in the right panel makes 
this less obvious. Our spectroscopic measurements, dis- 
played by the diamonds in Figure [3l have similar scatter 
to the measurements of MR95 and help to corroborate 
the presence of break, especially when using the N2 di- 
agnostic. The statistical analysis of the data of MR95 
suggests that for both metallicity diagnostics a model 
with a break provides a significantly better characteri- 
zation of the data. The same statistical tests performed 
on our spectroscopic data suggest that a model with a 
break provides only a marginally better fit to the data. 
The significantly lower confidence results from the sub- 
stantially fewer degrees of freedom in our fit owing to the 
greatly reduced sample size. 

The location of the break, in units of Rejj, is 1.0 and 
1.2 for the N2 and 03N2 diagnostics, respectively. The 
location of the break is in rough agreement between the 
two diagnostics and the value derived from the N2 diag- 
nostic is exactly consistent with MR95. For both metal- 
licity diagnostics, the slope of the outer gradient, within 
the errors, is consistent with zero. We do not expect 
the effective radius to play any significant role in deter- 
mining the gas dynamics driving the observed break in 
the abunda nce gradient. From numerical simulations of 
NGC 3359, iRozasI (|2008[ ) find that the corotation radius 
is 95". The effective radius for this galaxy is 87" and 
we attribute the break location of R^jj = 1 to the near 
coincidental location of the effective and corotation radii. 

5. IONIZATION PARAMETER AND THE R23 DIAGNOSTIC 

We determine the oxygen abundance and ionization 
parameter from the i?23 and 032 line ratios, respectively. 



032 = log{[OIII]/[OII]}. (10) 

Where [Oil I] is the sum of the hne flux of [0///]A4959 
and [O///]A5007 and [Oil] is the flux of the dou- 
blet centered at 3727A. This diagnostic builds on the 
photoion iz ation grids a nd stellar evolution models of 
iKewlev fcPopital (|200l . 

The i?23 method for determining meta llicities is known 
to depend on the ionization parameter. iPilvuginI (|2003[ ) 
has argued that artificial breaks in the abundance gra- 
dient may result from radial variations in the ionization 
parameter. The ionization parameter, which we denote 
by q, characterizes the ionization state of the gas within 
the HII region and quantitatively represents the number 
of ionizing photons per second per unit area divided by 
the hydrogen number density. The ionization parameter 
has the units of velocity and physically represents the 
maximum velocity of the ionization front through the 
nebula. The metallicity and ionization parameter are 
interdependent. We use equations A4 and A6 given in 
the appendix of IKewlev fc Ellisoiil ([2008) to determine 
the ionization parameter and metallicity respectively. In 
order to obtain a consistent measurement of the metal- 
licity and ionization parameter an iterative scheme is 
us ed, the details of which are provided in the appendix 
of IKewlev fc EllisonI (l200l . 

Metallicity is not a singly- valued function of i?23, but 
rather takes on two values for a given ratio. The peak of 
the R23 ratio occurs at 12 -f- log{0/H) - 8.4. This de- 
generacy is due to the fact that on the lower branch i?23 
increases with metallicity because the intensity of the 
coUisionally excited [Oil] and [OIII] lines increase with 
the abundance. On the upper branch, nebular cooling, 
which results from a coUisional excitation followed by 
photon emission, effectively cools the nebula decreasing 
the electron temperature leading to a decrease in the rate 
of coUisional excitation of the [Oil] and [OIII] lines. In 
order to break this degeneracy we use the [Nil] /[Oil] 
line ratio. We place objects on the lower branch of the 
R 23 diagnostic if loq([NII ]/[OII]) < -1.2 (see appendix 
of IKewlev fc EllisonI [2?)08h . 

Figure [4] shows the oxygen abundance and ioniza- 
tion parameter deterrnined f rom the iterative scheme of 
IKobulnickv fc KewlevI (f200l using the R23 and 032 line 
ratios, respectively, plotted against effective radius. We 
were able to determine the metallicities in only 13 of 19 
HII regions due to the non-detection of the [Oil] doublet, 
3 of which are found to be on the lower branch. MR95 
do not observe the [OH] line so we are unable to perform 
this diagnostic on their data. 

The abundance gradient shown in the top panel of Fig- 
ure [4] is suggestive of a break around R/Reff ^ 1. More 
importantly, there is no observed trend in the ionization 
parameter with radial distance. This provides strong ev- 
idence against the assertion that the break in the abun- 
dance gradient may be attributable to the radial varia- 
tions in the ionization parameter. 
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6. DISCUSSION AND CONCLUSIONS 

We build on the work of IMartin fc Rovl (jl995f) in order 
to validate the radial oxygen abundance gradient break 
they report. We find that the photometric determination 
of line ratios agrees reasonably well with the spectro- 
scopic line measurements. We make a small correction 
for systematic differences between the two methods of 
measurement which we attribute to uncertainties in con- 
tinuum subtraction and reddening corrections. Both the 
corrected data of MR95 and our spectroscopic measure- 
ments of 19 HII regions show evidence of a break in the 
gradient. We statistically compare a model with a single 
line to a model with a break. We find that the model 
with a break provides a significantly better characteriza- 
tion of the data. We determine the ionization parameter 
from the 032 line ratio and find that it appears to be 
u ncorrelated to th e radial distance. 

IPilvueinl ([200l examines the observed bends in the 
slopes of the radial oxygen abundance gradients in spiral 
galaxies. He argues that the upper/lower branch degen- 
eracy in the metallicity derived from i?23 along with gra- 
dients in the ionization parameter may result in a false 
break or bend in the radial oxygen abundance gradient. 
The main focus of this work is to address these issues. 

In order to address the i?23 degeneracy, we use an 
abundance indicator that is calibrated for these type 
of data, namely the spectroscopically determined N2 
line ratio measurements of individual HII regions in 
the metaUicity range 8.6 > 12 + log{0/H) > 8.2 
(jPettini fc Pagelll2004[ ). Because the N2 ratio IS a mono- 
tonic function of metallicity, it avoids the degeneracy is- 
sues that arise when using i?23. Additionally, for 13 of 
our 19 HII regions we were able to break the degeneracy 
in the i?23 diagnostic using the [NII]/[OII] line ratio. 
The abundance gradient determined from the i?23 diag- 
nostic appears to have a break consistent with the deter- 
mination from the N2 and 03N2 diagnostics, though it 
cannot be established from these data alone due to the 
small number of regions observed. 

The second issue raised by iPilvuginl (|2003l ) is that 
a radial oxygen abundance gradient break can result 
from a radial variation in the ionization parameter. He 
shows that for NGC 2403 a false bend results from 
using the R23 method for metallicity determination. 
We compile data from the literature for NGC 2403 
and determine the metallicity using the N2 diagnos- 
tic dMcCall et al.lll985l: iFierro et al.lll986t iGarnett et al.l 
119971: Ivan Zee et all I1998D . We find no evidence of a 
break. This suggests that the N2 diagnostic is less sus- 
ceptible to such confusion. Moreover, we examine the 
ionization parameter in NGC 3359 and find no appar- 
ent correlation with radial distance. We conclude that 
the break in the radial abundance gradient we observe in 
NGC 3359 is not due to a radial variations in the ioniza- 
tion parameter. In this way, we are able to statistically 
establish the abundance gradient break in NGC 3359 
while avoiding the problems associa ted with the points 
of concern raised bv IPilvuginl (|2003| ). 

The radially differential star formation rate is taken to 



be the cause of abundance gradients in spiral galaxies. 
However, it has been observed that barred galaxies gen- 
erally have shallower gradients due to large scale gas mix- 
ing associated with the bar formation. The presence of 
breaks in the abundance gradients is consistent with nu- 
merical simulations of barred galaxies with a young, s till- 
forming bar (Fried U et al.llT99llFViedh fc BendlT995[) . It 
is thought that because the bar is still forming and there 
is ongoing star formation, turbulent gas mixing has not 
had the time to smooth out the radial abundance gradi- 
ent inside the corotation radius, while outside the coro- 
tation radius, large-scale gas motions have homogenized 
the abundance. 

The formation and presence of a bar has been associ- 
ated w ith morphological transformations. iDutil fc Rovl 
(|1999f l argue that bars may be able to transform late- 
type galaxies into early-type ones. They find that the 
central abundance gradient of galaxies when plotted 
against morphological type is suggestive of two distinct 
sequences, one formed by late-type barred and early- 
type (barred and unbarred) galaxies and another com- 
prised of normal late-type galaxies. One interpretation 
of this sequence is that the bar flattens the abundance 
gradient and dissolves with time in early-type galax- 
ies. In simulations, the dissolution of the bar has been 
linked to the growth of the bulge and co uld be respon- 
sible for morphological tr ansformatio ns ([Combes et al.l 
119901 : iNorman etall 119961 |pfcnnigcr| \19m . However, 
the observational evidence for bar dissolution is limited 
(Bloc k et al . 2002; Da s et al. 2003). 

From recent infrared studies, the fraction of local spiral 
galaxies having bars is well esta blished and f ound to be 
around ^ 6 0% (Eskrid ge et al.l 120 00: Mar ino va fc Jogeel 
12007!; iMenendez-Delmestre et al.l i2007i) . iSheth et al. 
(2008) find that the bar fraction for spiral galaxies is 
greater for higher stellar mass systems and that this re- 
lation evolves from 0.2 < z < 0.8, with the greatest 
evolution in the fraction occurring for the lowest mass 
systems. This suggests that galaxy downsizing is not 
only a process of star formation regulation in massive 
systems but also one of dynamical evolution. 

The presence of bars in local spiral galaxies is preva- 
lent while the epoch of bar formation is short. Therefore, 
observations of galaxies where such abundance gradient 
breaks are expected to occur are rare. Bar formation 
may be a process leading to morphological transforma- 
tions and has been shown to exhibit downsizing. Abun- 
dance gradient breaks have been taken as evidence of the 
young bars still in the process of formation. Observing 
and understanding this transient epoch of formation may 
prove to be a very important aspect of understanding the 
processes governing the evolution of galaxies. 
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